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ABSTRACT 

Scaling relations that link asteroseismic quantities to global stellar properties are important for gaining un- 
derstanding of the intricate physics that underpins stellar pulsation. The common notion that all stars in an 
open cluster have essentially the same distance, age, and initial composition, implies that the stellar parameters 
can be measured to much higher precision than what is usually achievable for single stars. This makes clus- 
ters ideal for exploring the relation between the mode amplitude of solar-like oscillations and the global stellar 
properties. We have analyzed data obtained with NASA's Kepler space telescope to study solar-like oscillations 
in 100 red giant stars located in either of the three open clusters, NGC 6791, NGC 6819, and NGC 681 1. By 
fitting the measured amplitudes to predictions from simple scaling relations that depend on luminosity, mass, 
and effective temperature, we find that the data cannot be described by any power of the luminosity-to-mass 
ratio as previously assumed. As a result we provide a new improved empirical relation which treats luminosity 
and mass separately. This relation turns out to also work remarkably well for main-sequence and subgiant 
stars. In addition, the measured amplitudes reveal the potential presence of a number of previously unknown 
unresolved binaries in the red clump in NGC 6791 and NGC 6819, pointing to an interesting new application 
for asteroseismology as a probe into the formation history of open clusters. 
Subject headings: binaries: general — open clusters and associations: individual (NGC 6791, NGC 6819, 

NGC 6811) — stars: fundamental parameters — stars: interiors — stars: oscillations — 

techniques: photometric 
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1. INTRODUCTION 

The highly complex processes involved in the excitation 
and damping of stochastically excited (solar-like) oscillations 
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make estimation of their amplit udes from pulsa t ion modelling 
partic ularly challenging (e.g. iHoudekl [2006t ISamadi et all 
2007). A scaling relation for the amplitude has therefore 
been of sig nificant interest since it was first introduced by 
iKieldsen & Bedding! d!995l). Their 'L/M relation', based 
on the oretical work by Christensen-Dalsgaard & Frandsenl 
(1983) of near main-sequence stellar models, suggested that 
the amplitude in radial velocity would simply scale as the 
luminosity-to-mass ratio. Using observations of stars made in 
both radial velocity and intensity Kjeldsen & Bedding ( 1995) 
also suggested that the amplitude in intensity, A\, would scale 
as 

((L/L Q )/(M/M & )Y 

A " A/500nm(T off /5777K)'' A50tom -©' (i) 

where s = 1, A is the central wavelength of the photometric 
bandpass, and Asounm,© the observed solar value at 500nm. 
They found empirically r = 2.0, which was a slight modi- 
fication to r = 1.5 derived if they assumed the stellar oscil- 
lations to be purely adiabatic. Subs equen t modelling by e.g. 
iHoudeket all d 19991) : ISamadi et all d2007l) has lead to varia- 
tions of the L/M relation where, in essence, different powers 
of the L/M ratio h ave been derived (s = 0.7-1.3). Recently, 
IVerner et all (1201 ll) found s = 0.4-1.0 depending on T cS of 
a large sample (642) of main- sequence and subgiant stars ob- 
served by Kepler dKoch et alj|2010h . 

The existence of solar-like oscillations in red giant 
stars is now well establi shed ob servationally, most re- 
cently fr om CoRoT (e.g. [d e Ridder et al. 2009J) and Ke- 
pler (e.g. Gilliland et al. 2010; Bedding et al. 2010b), as well 
as theoretically dDupret e t al. 2009; Montalban et all 120101 
iDi Mauro et al.1 1201 lb . Despite the significantly different 
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structures of red giants compared to the stars and models 
on which the L/M relation has been founded, the absence 
of an alternative has also seen this relation widely used for 
red giants, including s everal attempts to determine the best 
matching exponent , s (IStello et alJl2007t lMosseretal"1l2010t 
Baudi net al.ll201 lh . While the majority of results on red gi- 
ants are on field stars, the recent clear detections in open clus- 
ter red giants emerging from Kepler dStello et al.l l2010.Paper 
I) has opened up the seismic exploration of clusters and the 
advances that clusters bring to the interpretation of asteroseis- 
mic data ( Basu et al.l201 luHekker e t al. 201 1; Miglio et al. in 
prep.; Stello et al. submitted). In particular, stars in an open 
cluster are thought to share a common distance and initial 
chemical composition, which allows one to derive the stellar 
luminosity to much higher precision than for most field stars. 
In addition, the common age of red giant stars within each 
cluster implies that they have practically the same mass, re- 
sulting in a relatively low uncertainty on their measured mean 
mass assuming there is no significant mass loss (Miglio et al. 
in prep.). Combined with high quality standard photometry 
we can therefore obtain more robust predictions of the am- 
plitudes from scaling relations and hence investigate these in 
ways not possible for the field stars observed by the current 
space mission CoRoT and Kepler. 

Based on only one month of Kepler data of a single cluster, 
in Paper I we already demonstrated the potential for inves- 
tigating the L/M relation by taking advantage of the com- 
mon cluster properties of the stellar sample. We now have 
Kepler time-series photometry that span 10 times longer for 
stars in three open clusters (NGC 6791, NGC 6819, and 
NGC 6811), which exhibit distinctly different stellar masses. 
In this paper we are therefore extending considerably the anal- 
ysis of the amplitude scaling relation for solar-like oscilla- 
tions. 

2. OBSERVATIONS, TARGET SELECTION & CLUSTER 
PARAMETERS 

The photometric time-series data were obtained between 
2009 May 12 and 2010 March 20 (observing quarters 1-4), 
providing approximately 14,000 data points per star obtained 
in the spacecraft's long-cadence mode (Ai = 29.4 min). A 
detailed description of the data reduction from raw images to 
final li ght curves is given in Jenkins et al.l (l2010l) : lGarcia et all 
(12011b and Stello et al. (submitted). 

Our initial star sample was the one selected by Stello et 
al. (submitted), who used seismic and conventional mea- 
surements to identify cluster membership and blending of 
each star. We excluded the seismic non-members and further 
trimmed the sample by removing the brightest (largest) and 
faintest stars, for which the measurement of the mode ampli- 
tude would not be reliable due to: (1) difficulty in determin- 
ing the noise level at low frequency in the power spectrum of 
the largest stars (oscillating at very low frequencies) and (2) 
low signal-to-noise and potential blending of the faintest stars. 
The increased flux in the photometric aperture from a blend- 
ing star, such as an unresolved binary companion, will tend 
to reduce the relative flux variation that we measure as the 
oscillation amplitude. To minimize this bia s further, we ex- 
cluded a total of 23 spectroscopic binaries dHole et al.ll2009l) 
and stars that we expected to be binaries based on their lo- 
cation in the color-magnitude diagram. This still left a large 
sample of 100 stars for further analysis. We finally investi- 
gated effects of blending of single stars based on the results 
by Stello et al.. Few of the blended stars indicated by Stello et 
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FIG. 1. — H-R diagram of the selected cluster stars. Small gray sym- 
bols mark the known and potential binaries. Representative isochrones from 
IMarigo et alj [2008) (NGC 6791: 5.6 Gyr Z = 0.030 and NGC 6819: 
2.4 Gyr, Z = 0.019) and Pietrinfer nTeTaT] J20041) (NGC 6811: 500 Myr, 
Z = 0.008) are shown to guide the eye. 



al. showed lower than expected amplitudes, but no rigorous 
criterion for when blending had a significant impact on the 
amplitude could be obtained from those results. We therefore 
did not exclude any of our remaining stars that were listed as 
blends. 

We adopted the luminosities, masses and effecti ve temper- 
atures from Stello et al. ( submitted). We refer to iBasu et al.1 
d2011l) and iHekker et ail d2011l) for further details on the 
derivation on the mass and effective temperature, respectively. 
In summary, the average mass (here adopted for each star) is 
1.20 ± 0.01 M (NGC 6791), 1.68 ± 0.03 M Q (NGC 6819), 
and 2.35 ± 0.04 M (NGC 681 1), while the luminosities and 
temperatures of our final sample are shown in Figure Q] and 
have typical uncertainties of ~ 10% and ~ 2%, respectively. 

3. MEASUREMENT OF OSCILLATION AMPLITUDES AND i/ max 

Oscillation amplitudes were extra cted by five d ifferen t 
teams using pipelines descr ibed in iHekkeretalJ (1 2010): 



Hube retalJ (120091): iKallinger et al.1 d2010l): iMathur et alj 
(2010): lMosser & Appourchaux (2009). These methods are 
all based on the measurement of the integrated oscillation 
power, which we converted to an amplitude per radial mode. 
The integrated power was fou nd either by smoothi ng the 
power spectrum as described by Kield sen et alj d2008l) or by 
fitting a Gaussian function to the oscillation power envelope. 
Figure|2]shows the former. To obtain the amplitude per radial 
mode the oscillation power (P bs, Figure |2]l is multiplied by 
Av to obtain the power per radial order, where Av is the fre- 
quency separation between consecutive radial orders. Finally 
we divided by the f actor, c, wh ich is the effective number of 
modes per Av (Kieldsen et al. 2008). We adopted the so- 
lar value c = 3.04 from Bedding et al. (2010a), which agrees 
well with the measured mean value for red giants (Mosser et 
al. submitted). We note that our final results (Sect. [4731 were 
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FIG. 2. — Power spectrum of a typical star. The smoothed spectrum (solid 
white line) and fit to the stellar granulation background (dashed line) are 
shown. The oscillation power, P b s > is evaluated at the frequency of max- 
imum power, c max . 



not affected signif icantly if we adopted the recent factor by 
Ballot etail d2011l) . Hence, the observed amplitude per radial 
mode, A \> s (l = 0) was derived as: 



A ohs (l = 0) = (PobsAi//3.04) 



1/2 



(2) 



For this we normalized the power spectra according 
to the amplitude-scaled v ersion of Parseval's theorem 
(Kjeldsen & Frandsen 1992), in which a sine wave of ampli- 
tude, A, provides a peak in the power spectrum of A 2 . The 
typical uncertainty in the measured amplitude is ~ 10%. 

To explore whether the applied solar conversion factor, c, 
provided reasonable amplitudes for red giants, we ran sim- 
ulations that as input too k pulsation frequencies deriv ed us- 
ing the ADIPLS code (Christensen-Dalsgaard 2008a) for a 
representative set of ASTEC models (Christensen-Dalsgaard 
2008b). Details of th e simulator can be found inlChaplin et al.l 
(2008). Following Christensen Dalsgaardl ( 120041) . the input 
mode amplitudes were scaled relative to the radial modes us- 
ing the mode inertia, /, as A oc I~ 2 . Despite significant 
differences in the frequency spectra of red giants compared 
to the Sun, in particular the presence of many mixed mode s 
(iDupret et all l2009t iBlicketail l201lt iBedding et al.l l201lh . 
the pipelines returned amplitudes within 10% of the input val- 
ues. We regard this as acceptable given the uncertainty from 
intrinsic scatter of the oscillations and the slightly different 
approaches for extracting the amplitudes in each pipeline, in 
particular the fitting and subtraction of the stellar granulation 
background (Mathur et al. submitted). Based on a representa- 
tive set of stars, we found good agreement between the differ- 
ent pipelin es. In this p aper we show the results from the S YD 
pipeline dHuber et al.l 120091) . which provided amplitudes for 
the widest range o f stars, and we compar e our final result with 
the CAN pipeline dKallinger et alj|2010b . which exhibited the 
largest overlap in stellar sample with the SYD pipeline. Both 
pipelines show robust performances in their estimation of the 
stellar granulation backgro und (Mathur et al. submitted). We 
refer to lVerneret al.l (12011 ) and Mosser et al. (submitted) for 
detailed amplitude comparisons. 

In addition to amplitude, the pipelines also measured the 
frequency of maximum power, u max (Figure |2). The uncer- 
tainties in ^ max are typically 1-2%. 

4. RESULTS 
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Fig. 3. — Observed amplitude versus f max for stars in NGC 6791 
(red diamonds), NGC 6819 (purple triangles), and NGC 6811 (blue 
squares). The binary stars are shown with small gray symbols. The 
clump stars are marked. The dashed line shows a power law with slope 
—0.75. Colored lines are the cluster isochrones (Figure [T} where ampli- 
tude and lAnax have been derived using Equation [T] with s = 0.75 and 
i/ max = (M/M Q )/(L/L Q )(r cff /5777K)3- 5 3100MHz. The black cross 
at (10,100) indicates a typical 1-cr error bar. 



4.1. A, 



bs versus v n 



As noted by IStello et alj (12007b : iMosser et alj ( 120101) : 
iHuber et all d201 0). it can be convenient to plot the measured 
amplitude as a function of v max , since the currently adopted 
scaling relations predict a simple relation between the two. In 
particular, by dividing A\ oc {L/M) s T~ ff r (Equation [B by 
(i/ mai ) s ex (M/L) a T*g 5a dBrownetal.1119911) and rearrang- 



ing, we obtain A\ oc i/ m * x T ( 



3.5s- 



Hence, such a purely em- 



pirical plot allows one to make some inference on how the am- 
plitude depends on the stellar parameters L, M, and T e g even 
when t hose are not very well known (e.g. IMosser et alj|20ldl : 
Huber et al. 2010; Huber et al. submitted; Mosser et al. sub- 
mitted). 

In Figure [3] we show the measured amplitude as a func- 
tion of f max , where each set of symbols present results of one 
cluster. We also mark the location of the clump of helium- 
core burning stars for each cluster, which illustrates the large 
range in v max arising mainly from the difference in the stellar 
mass between the clusters. 

Guided by the fiducial dashed line, we see that stars within 
each cluster roughly follow a power law with exponent —0.75, 
but with a clear offset from one cluster to another by up to 
~ 50%. The more massive the stars, the lower the oscilla- 
tion amplitudes at a given v max . This offset is not expected 
from the scaling relations for A\ and ^ max , as illustrated by 
the isochrones in Figure [3] Since the scaling relat ion for 
i/ mal is probably good to within a few percent ( Stell o et all 
2009; Bel kacem et alJ 1201 ll) . the observed offsets strongly 
suggest that (L/M) 8 T~ s r does not adequately predict the am- 
plitude for these sta rs. From a large sample of field red giants 
IHuber et all d2010h noted that the scatter in the amplitude at 
a given i/ max was larger than expected from the uncertain- 
ties and that this indicated a spread in mass in their sample. 
However, a qualitative analysis was not attempted due to the 
relatively large uncertainties in the fundamental stellar param- 
eters. Fortunately, with our cluster sample we can directly 
fit the measured amplitudes to their predictions derived from 
well-constrained stellar parameters. 
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4.2. Fitting the L/M relation 

First, we fitted the observed amplitudes to the predicted 
amplitudes for NGC 6819. For this purpose we derived the 
predicted amplitude using the L/M relation (Equation Q~|i 
and adopting A = 650 nm as the central wavelength of 
the Kepler bandpass, hence ^4 b s .© = 3.49 (peak scaled) 
dMichel et a l. 2009). The least-squares fit resulted in s = 
0.76±0.01 when adopting the empirical value of r = 2, which 
is the value of r we will adopt in the following. Using r = 1.5 
only has the effect of increasing s by about 0.03. This result 
is compatible with Paper I, which qualitatively found the best 
match for s to be slightly higher than 0.7. When repeated for 
NGC 6791, we found s = 0.87 ± 0.01. The small number of 
stars in NGC 6811 did not merit a fit on its own, but the two 
other clusters already indicate inconsistent results. 

Hence, we tried next fitting all three clusters simultane- 
ously. Due to the correlation between M and T e g (the hotter 
and younger clusters have more massive stars; Figure [TJ, we 
still kept r fixed. Figure HJa) shows the result. The best fit 
resulted in s = 0.74 ± 0.01. It is apparent that the clusters 
are offset from one another, as expected from Figure [3] but 
we also see that the fit systematically underestimates the am- 
plitude for the most luminous stars. If r was treated as a free 
parameter we did obtain a better fit overall, but it still un- 
derestimated the amplitudes of the stars in NGC 6791, and 
in particular the most luminous stars in the sample, by 20- 
30%. In summary, while the (L/M) s scaling provided ac- 
ceptable results when fitted to one cluster at a time (although 
giving different results for s), our analysis has demonstrated 
that (L/M) s cannot explain the observations in all clusters 
simultaneously. 

4.3. A new scaling relation for amplitudes 

In the following, we therefore fitted the exponents on L 
and M independently, hence A\ oc L s M~ t T~^' . The re- 
sult, shown in Figure E|b), is a much improved fit where all 
three clusters fall on top of each other and follow the one-to- 
one relation. The best fitting parameters are s = 0.90 ± 0.02 
and t = 1.7 ± 0.1 - the same as we obtained from first con- 
verting A bs to a bolometric amplitude (Ballot et al., submit- 
ted) and then fitting to A bo \ oc L S M~*T^. For the stars 
with ,4 obs £ 80ppm the scatter of A obs /(L - 90 Af- 1 - 7 T cff 2 ) 
is 14%, in perfect agreement with the quoted uncertainties on 
A bs, L, M, and T c g. The increased scatter (22%) towards 
lower luminosity stars is potentially due to remaining issues 
of blending in the sample and/or an increase in the uncertain- 
ties of the measured amplitudes for the faintest stars. The 
latter was, however, not reflected in the estimated uncertain- 
ties reported by the pipelines showing only slightly increased 
uncertainties at most. Again, under the adiabatic assumption 
(r = 1.5) s would slightly increase (to 0.95) as would t (to 
1.8). 

To investigate the robustness of our fit we did the following. 
If we ignored the NGC 6811 stars in the fitting, the result and 
hence the excellent alignment of all three clusters was very 
similar (s and t within la). This is perhaps not surprising 
given the few stars in our NGC 6811 sample. Nevertheless, 
this result is reassuring since the amplitudes of NGC 681 1 are 
then correctly predicted from a fit based only on NGC 6791 
and NGC 6819. We further investigated the effect on the fit 
if we ignored all clump stars to obtain an even more homo- 
geneous sample, which showed practically no change to the 
best fitting parameters. This indicates that any possible mass 
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FIG. 4. — (a) Observed versus predicted amplitude for the best fitting rela- 
tion of the form A\ oc (L/M) S T~ S . Symbols are the same as in FigurefS] 
Binaries, which are shown with small gray symbols, were not included in the 
fit. (b) As panel (a) but fitting to A x oc L S M ~ l T~g . The inset shows the 
X 2 near its minimum, (c) Illustration of how well the fit in panel (b) predicts 
amplitudes for other main-sequence, subgiant, and red giant stars (see text). 
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loss, which is expected to occur predominantly near the tip 
of the red giant branch, has no effect on our result. A small 
systematic change of a few percent on s and t was, however, 
observed by removing some of the most deviant stars at low 
amplitudes. Finally, we repeated the fit on the sample of stars 
that were in common between the SYD and CAN pipelines. 
The differences in s and t based on these different pipelines 
were 2% and 15% in s and t, respectively, the latter only just 
within 3<r of the formal uncertainty. 

We finally tes ted the new scaling relation suggested by 
Kjeldsen & Bedding (2011), but found it to overestimate the 
amplitude for the cluster stars similar to the result found by 
Huber et al. (submitted) and Mosser et al. (submitted). 

4.4. Main-sequence and subgiant stars 

Now, with an improved scaling relation for red giant stars, it 
is interesting to see how well it applies to main-sequence and 
subgiant stars. To investigate this we took amplitude measure- 
ments of the Kepler field stars presented by Huber et al. (sub- 
mitte d), the CoRoT F-ty pe stars HD49933 and HD 18 1420 
from lMichel et al] (I2008J) (conv erte d to A n -h R (l=Q)\ an d Pr °- 



cyon from Arentoft etafl d2008l) and lHuberet ail (120111) . The 
amplitude measurement in velocity of Procyon was converted 
to intensity using models by Houdek (2010). We used our new 
scaling relation to predict the amplitudes based on L, M, and 
T e g fr om Huber et al. (submitt ed) (Kepler sample), Bruntt 
d2009h (HD49933/181420), and iBonanno etail (120071) (Tro- 
cyon). Given that the new relation is only based on the cluster 
red giants, it is remarkable how well it agrees for this broad 
range of stars (Figure[4j;). We note that the uncertainty in the 
mass of the Kepler (~ 10-20%) and CoRoT (~ 5-10%) field 
stars is significantly larger than for Procyon (~ 2%) and the 
cluster stars (~ 1-2%). While the values of s and t found in 
this Letter are slightly different from, although still in agree- 
ment within the uncertainties, those found by Huber et al. 
(submitted) for the Kepler field stars, the qualitative agree- 
ment across all stars is quite similar to that found by Huber et 
al. (their Fig. 5). 

4.5. Unresolved binaries 

It is evident, particularly from Figure|4|b), that many of the 
known and potential binaries (small gray diamonds and trian- 
gles) show relatively low amplitudes. For NGC 6791 we had 
no spectroscopic determination of binaries, but a significant 
fraction of its red clump stars show lower than expected am- 
plitudes and hence strong evidence for 'diluted' light curves 
due to the presence of unresolved binary companions. This 
shows a new exciting way of applying asteroseismology to 



identify binary stars and hence to probe the formation of these 
stars in clusters, which will be investigated in detail in a forth- 
coming paper. 

5. CONCLUSIONS 

Our analysis of solar-like oscillations in 100 red giant stars 
in three open clusters revealed that previously adopted scaling 
relations based on the luminosity-to-mass ratio for predicting 
amplitudes are not adequate for red giants. We found an em- 
pirical scaling relation by fitting the observed amplitudes to 
a more general form than the previous L/M relation. The 
result, 



AxKlP-w/iM 1 - 7 ^) 



(3) 



and 



4,oi«L a90 /(M 1 - 7 T eff ), (4) 

which showed considerable improvement for red giants, 
turned out to also work remarkably well for main-sequence 
and subgiant stars. 

Interestingly, the lower than expected amplitudes of some 
red clump stars in NGC 6791 and NGC 6819 revealed that 
they were likely unresolved binaries, many of which were not 
known previously. This method for identifying binaries could 
add interesting new insight to the formation history of these 
clusters. 

In this investigation we ignored an y possible effect on a m- 
plitude from metallicity differences (Sama di et alJ 120071) of 
the three clusters, which have value s of [Fe/HlNOCfi 73i — 0.3 
and [Fe/H]NGC68i9 — 0.1 (see iBasu et all 1201 ll and refer- 
ence therein), while it is unknown for NGC 6811. To improve 
on that will require better determination of the cluster metal- 
licities. In addition, we would need more clusters (with sig- 
nificantly different stellar parameters) to allow the fitting of a 
rigorous empirical relation including one more free parameter 
such as metallicity. 

With more Kepler data in the future, we expect to have clus- 
ter stars covering a large range of T e g, which will include 
turn-off stars at the end of the main sequence, allowing us to 
also fit the exponent, r, of the T c g dependence in the ampli- 
tude scaling relation. 
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